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摘要 : 全 球 卫 星 导 航 系 统 空间 信号 分 析 是 系统 设计 阶段 与 运行 过 程 中 的 一 项 重要 工作 。 
首先 给 出 了 基于 全 球 卫 星 导航 系统 监测 接收 机 观测 量 的 空间 信号 稳定 性 分 析 方 法 ， 然 后 通过 
数字 仿真 实验 证 实 了 方法 的 可 用 性 ， 最 后 以 北斗 系统 地 球 同 步 轨道 卫星 的 实测 伪 距 数据 为 
例 ， 从 原始 伪 距 观测 数据 、 伪 距 拟 合 残 差 与 伪 距 拟 合 残 差 标准 差 3 方面 对 信号 稳定 性 进行 了 
分 析 ， 实 验 结果 进一步 验证 了 本 文 方法 的 可 行 性 与 有 效 性 。 方 法 对 于 保障 卫星 导航 信号 的 连 
续 性 与 可 靠 性 具有 一 定 意义 。 
关键 词 ， 全球 卫星 导航 系统 ; 信号 质量 分 析 ; 稳定 性 ; 伪 距 
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随 着 卫星 导航 技术 的 发 展 ， 世 界 各 大 国 相继 建设 了 各 自 的 全 球 卫 星 导 航 系统 ( Global Navigation 
Satellite System ,GNSS) 或 区 域 卫星 导航 系统 ， 例 如 美国 全 球 定位 系统 、 俄 罗斯 GLONASS 、 欧 盟 伽 利 
o 略 、 中 国 北斗 系统 。 全 球 卫星 导航 系统 的 空间 信和 号 质量 不 仅 直接 影响 着 卫星 导航 系统 定位 、 导 航 和 授 
= 时 ( Positioning, Navigation and Timing, PNT) 的 性 能 ， 而 且 能 够 反映 卫星 有 效 载 符 的 在 轨 状 态 与 各 种 电 
^ 性 能 指标 ， 因 此 ， 全 球 卫 星 导航 系统 空间 信号 质量 是 衡量 卫星 导航 系统 优 劣 的 一 个 重要 指标 ， 相 应 的 
- 信和 号 质量 评 佑 技术 可 以 为 信号 设计 、 信 和 号 在 轨 测 试 以 及 信和 号 在 轨 监 测 提供 重要 支撑 ， 从 而 为 系统 故障 
诊断 和 空间 信号 完好 性 监测 等 提供 重要 保障 … 。 

由 于 全 球 卫星 导航 系统 空间 信号 评估 是 系统 设计 阶段 与 运行 过 程 中 的 一 项 重要 工作 ， 因 而 引起 了 
各 国 的 高 度 重 视 ， 也 是 当前 卫星 导航 领域 的 一 个 研究 热点 。 目 前 ， 全 球 定位 系统 、 北 斗 系统 和 伽利略 
等 卫星 导航 系统 均 建 立 了 配套 的 空间 信号 质量 监测 和 分 析 系 统 ， 如 全 球 定位 系统 由 美国 斯 坦 福 国际 研 
究 院 建立 了 导航 信号 监测 和 特性 分 析 系 统 ， 该 系统 利用 47 m 大 口径 抛物 面 天 线 对 全 球 定位 系统 空间 
信和 号 进行 全 天 候 实时 监测 扣 ;伽利略 由 欧洲 空间 技术 研究 中 心 导 航 实验 室 建设 了 导航 信号 特性 分 析 
系统 ， 该 系统 由 全 向 天 线 、 标 准 测量 仪器 和 测试 接收 机 等 设备 组 成 ”; 北斗 系统 由 中 国 科学 院 国家 
授时 中 心 搭建 了 全 球 卫星 导航 系统 空间 信号 分 析 平 台 ， 以 标准 测量 仪器 分 析 、 离 线 数 据 评 估 与 监测 接 
收 机 3 种 方式 对 包括 北斗 系统 空间 信和 号 在 内 的 多 种 全 球 卫星 导航 系统 信号 质量 进行 监测 和 评估 …” 。 

离线 数据 分 析 一 方面 根据 基带 信号 分 析 卫 星 在 短 时 期 内 的 工作 状态 5 ， 另 一 方面 通过 监测 接收 
机 测量 数据 对 卫星 的 长 期 工作 性 能 进行 评估 '“”。 本 文 主要 研究 基于 监测 接收 机 观测 量 的 卫星 导航 信 
号 稳定 性 分 析 方 法 。 


1 卫星 导航 信号 构成 
通过 简要 描述 全 球 卫星 导航 系统 发 射 的 信号 大 致 了 解 卫星 导航 信号 的 组 成 与 特点 。 导 航 卫星 发 射 


的 信号 从 结构 上 可 分 为 测 距 码 ( 伪 随 机 码 ) 、 载 波 与 数据 码 三 个 层次 所 。 在 卫星 发 射 信号 之 前 ， 先 对 数 
据 码 与 测 距 码 进行 调制 ， 然 后 再 将 调制 信号 和 正弦 波形 式 的 载波 混 频 ， 最 后 卫星 把 调制 后 的 载波 信号 
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发 射出 去 。 图 1 描述 了 导航 卫星 信号 的 基本 构成 。 
在 卫星 信号 接收 终端 ， 全 球 卫 星 导 航 系 统 接 
收 机 根据 接收 的 卫星 信号 输出 伪 距 (依据 测 距 码 
相位 测量 得 到 的 伪 距 也 称 作 测 码 伪 距 ) 与 载波 相 
位 (依据 载波 相位 测量 得 到 的 伪 距 也 称 作 测 相 伪 
WB) 两 个 基本 观测 量 。 下 面 介绍 伪 距 和 载波 相位 
两 个 基本 观测 量 的 观测 方程 。 
全 球 卫 星 导航 系统 接收 机 要 实现 定位 、 导 航 
或 者 授时 ， 必 须 准确 测量 卫星 与 接收 机 之 间 的 距 
离 。 伪 距 和 载波 相位 是 接收 机 对 观测 卫星 的 两 个 
基本 距离 观测 量 。 伪 距 和 载波 相位 一 旦 发 生 中 断 或 者 跳 变 ， 定 位 、 导 航 和 授时 精度 都 不 可 避免 地 受到 
影响 。 伪 距 P 和 载波 相位 工 的 观测 方程 可 分 别 表示 为 
Psp +c(6t, -8) +T +È +e p, (1) 
Lep +e(ôt, 00) +T —D AM +e. (2) 
T 其 中 , p: 为 卫星 信号 发 射 时 刻 的 卫星 位 置 到 信号 接收 时 刻 的 接收 机 位 置 的 距离 ，6t, 和 OD 分 别 为 接收 
2 机 时 钟 钟 差 和 卫星 时 钟 钟 差 ; T 和 天 分 别 为 对 流 层 折射 误差 和 电离 层 折 射 误 差 ; NI 为 载波 相位 的 整 
周 模糊 度 ; A 为 信号 波长 ; c 为 光 在 真空 中 的 传播 速度 ; si 为 伪 距 观测 噪声 和 多 路 径 效应 误差 之 和 ; 
a6 2 ,为 载波 相位 观测 噪声 和 多 路 径 效应 误差 之 和 。 


二 2 基于 全 球 卫星 导航 系统 接收 机 的 信号 稳定 性 分 析 方 法 


图 1 导航 卫星 信号 构成 示意 图 


Fig.1 Diagram of composition of navigation satellite signals 


is 全 球 卫星 导航 系统 空间 信号 的 稳定 性 分 析 一 般 通 过 评估 监测 接收 机 输出 的 伪 距 与 载波 相位 的 稳定 
c 性 实现 。 全 球 卫星 导航 系统 信号 稳定 性 分 为 短期 稳定 性 (1 天 或 几 小 时 ) 和 长 期 稳定 性 (一 星期 ) 。 信 
N 号 短期 稳定 性 主要 受 接 收 机 短暂 失 锁 和 接收 机 时 钟 钟 跳 等 因素 的 影响 ， 这 些 因素 导致 载波 相位 和 伪 距 
2 观测 值 突然 产生 跳 变 ， 然 后 又 恢复 正常 ， 而 信号 长 期 稳定 性 主要 与 卫星 时 钟 有 关 ， 如 卫星 钟 频率 漂移 
导致 载波 相位 和 伪 距 观测 量 出 现 缓慢 变化 ， 致 使 两 个 观测 量 整 体 产 生 漂移 ， 这 种 现象 反映 了 信号 的 长 
期 稳定 性 '"。 本 节 介绍 全 球 卫 星 导 航 系统 信号 稳定 性 的 分 析 方 法 ， 为 评估 载波 相位 和 伪 距 测量 值 的 
稳定 性 ， 主 要 分 析 二 者 的 曲线 拟 合 残 差 。 首 先 利用 一 段 时 长 的 载波 相位 或 伪 距 测量 值 ， 对 其 进行 多 项 
- 式 拟 合 ， 然 后 获得 多 项 式 拟 合 残 差 。 多 项 式 拟 合 依据 最 小 二 乘 原理 对 原始 观测 数据 进行 拟 合 ， 使 得 原 
一 台 观 测 数 据 和 拟 合 曲 线 之 间 具 有 最 小 的 差 值 。 具 体 实现 方法 为 ， 给 定 观 测 数据 (x,, y,) (i=1, 2,…， 
m), SHE n TER f (x), ARE ri my, f(x) (i81, 2, 5, m) 的 平方 和 最 小 ， 即 


m 


E- » [; * (2  - Y, P -y = min, (3) 


满足 (3) 式 的 (x)= Zar 称 作 最 小 二 乘 拟 合 多 项 式 ， 其 中 Go 1，2，…， n) 为 多 项 式 拟 合 参数 。 从 
几何 角度 而 言 ， 最 小 二 乘 拟 合 的 目的 是 寻找 与 给 定点 (*, y) (1= 1, 2, …, mn) 距离 平方 和 最 小 的 曲线 
y, f (x) (i81, 2, 5, m), Wf roy, - f i) (i=1, 2,…, m) 被 称 作 拟 合 残 差 。 对 于 载波 相位 和 伪 距 
拟 合 而 言 ， 电 离 层 折射 误差 和 对 流 层 折射 等 缓慢 变化 误差 已 经 大 体 消除 ， 因 此 拟 合 残 差 能 够 较 好 地 反 
应 多 项 式 拟 合 的 好 坏 。 实 际 上 ， 根 据 最 小 二 乘 原理 获得 的 最 终 拟 合 结果 可 以 充分 消除 观测 数据 中 的 系 
统 误差 ， 同 时 也 能 在 一 定 程度 上 削弱 随机 误差 的 影响 ， 因 此 得 到 的 拟 合 结果 具有 较 高 的 可 靠 性 。 
根据 多 项 式 拟 合 残 差 的 变化 情况 评估 载波 相位 和 伪 距 测量 值 的 稳定 性 ， 基 本 步骤 如 下 : 
(1) 读 取 全 球 卫 星 导航 系统 监测 接收 机 的 原始 观测 数据 ， 包 括 载波 相位 观测 值 和 伪 距 观测 值 ; 
(2) 绘制 载波 相位 观测 值 和 伪 距 观测 值 曲线 图 ， 根 据 曲 线 图 判断 观测 数据 中 是 否 存 在 跳 变 点 ; 
(3) 截 取 一 定时 长 的 连续 观测 数据 进行 多 项 式 拟 合 ， 获 得 多 项 式 拟 合 残 差 ; 
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(4) 根 据 多 项 式 拟 合 残 差 的 变化 范围 分 析 观 
测 值 的 稳定 情况 。 

图 2 为 分 析 载 波 相 位 和 伪 距 观测 量 稳 定性 的 
算法 流程 图 。 


3 ”数字 仿真 


3.1 方法 的 可 行 性 验证 

1 天 的 原始 伪 距 和 载波 相位 观测 数据 的 曲线 
图 近似 正弦 或 余弦 荫 数 (可 以 参见 第 4 节 给 出 的 
原始 伪 距 测量 曲线 ) 。 为 不 失 一 般 性 ， 首 先 产生 
一 个 周期 (0~2 普 ) 的 三 加 高 斯 白 噪声 的 余弦 函数 
以 仿真 伪 距 观测 值 ， 总 共 629 个 数据 点 ， 如 图 3， 
然后 分 别 截取 0~0. 6 m 5; 0.7 0-1. 3 m 这 两 段 模 
拟 数据 进行 拟 合 ， 这 两 段 数 据 均 包含 拐点 。 表 1 
列 出 了 释 加 不 同 高 斯 白 品 声 情况 下 的 多 项 式 拟 合 
残 差 标准 差 ， 记 为 o,， 同 时 表 1 也 给 出 了 未 县 
加 高 斯 白 噪声 时 的 拟 合 残 差 标准 差 ， 记 为 ci 。 
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14 d 


读 取 监 测 接收 机 原始 观测 量 
〈 伪 距 和 载波 相位 ) 


绘制 拟 合 残 差 曲线 ， 
拟 合 残 差 标准 差 


图 2 基于 全 球 卫星 导航 系统 接收 机 的 伪 距 和 
载波 相位 稳定 性 分 析 流 程 图 
Fig.2 Flowchart of the method for analyzing the stability 


of pseudo-range and carrier-phase 


R1 多 项 式 拟 合 残 差 标准 差 r 


Table 1 Standard deviation of polynomial — 0-0.6n 

T- 一 一 077 1-1.3m QS 
— fitting residuals Pa 
. d 
~ 高 斯 0, T, F / 

HRE —0-06 0.7-L3 | 0-06 0.7-1.3 9 t y 

标准 差 m " m m f 

^ ES ka 

0. 020 0. 001 0. 002 0.004 — 0.020 | E 2 

0. 040 0. 001 0. 002 0.044 0.040 | A 

0. 050 0. 001 0. 002 0.053 0.050 


从 表 1 可 以 看 到 ， 对 于 未 大 加 高 斯 白 噪声 的 
= 数据 拟 合 ， 多 项 式 拟 合 残 差 标准 差 o, 均 比 较 低 ， 
接近 于 0， 这 说 明 运 用 多 项 式 拟 合 方法 获得 的 拟 
合 结果 与 原始 观测 值 之 间 符 合 得 较 好 。 以 未 过 加 
高 斯 白 噪 声 的 多 项 式 拟 合 残 差 标准 差 为 参考 ， 向 余弦 函数 中 加 入 标准 差分 别 为 0.020 0.040 和 0. 050 
的 高 斯 白 噪 声 ， 获 得 了 含 品 数据 与 拟 合 数据 之 间 的 拟 合 残 差 标准 差 o,。 经 过 对 比 可 以 发 现 ，o, 和 高 
斯 白 噪 声 的 标准 差 非 常 相近 甚至 相等 ， 这 进一步 表明 ， 多 项 式 拟 合 方法 对 含 噪 数据 的 拟 合 结果 与 原始 
数据 之 间 也 具有 很 好 的 符合 性 ， 剩 余 的 拟 合 残 差 一 般 是 混 匠 在 原始 数据 中 的 随机 噪声 成 分 。 

为 更 加 直观 地 显示 多 项 式 拟 合 结果 ， 图 4(a) 和 图 4(b) 分 别 绘 出 了 0-0.6m 50.7 m-1.3 a 两 段 
含 噪 数据 曲线 和 多 项 式 拟 合 曲线 之 间 的 比较 。 从 图 4 可 以 看 出 ， 多 项 式 拟 合 方法 在 图 4(a) 和 图 4(b) 
曲线 拐点 处 的 拟 合 效 果 不 佳 ， 但 是 总 体 而 言 ， 图 4(a) 和 图 4(b) 曲线 的 整体 拟 合 残 差 较 小 ， 拟 合 残 差 
标准 差 与 琶 加 的 高 斯 白 噪声 的 标准 差 非 常 接近 ， 可 见 曲线 拐点 对 多 项 式 拟 合 效果 影响 不 大 ， 所 以 在 曲 
线 拐 点 处 ， 多 项 式 拟 合 仍然 有 效 。 

3.2 观测 异常 仿真 

利用 数字 仿真 分 析 观 测 数据 异常 条 件 下 的 多 项 式 拟 合 效果 ， 包 括 卫 星 时 钟 钟 跳 等 导致 的 测量 数据 

跳 变 、 接 收 机 失 锁 等 因素 造成 的 观测 数据 缺失 。 首 先 采 用 0~2 a 一 个 周期 的 含 噪 余弦 曲线 仿真 1 天 


图 3 dm EE E HS (HAZA 0.020) 的 余弦 曲线 
Fig.3  Cosine function with White Gaussian Noise 


(WGN) (00.020) 
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的 载波 相位 或 伪 距 观测 数据 ， 总 共 6 284 个 数据 点 。 然 后 分 别 在 0.13 0-0. 15 c BERI O. 76 o 点 模拟 接 
收 机 失 锁 等 引起 的 数据 缺失 和 卫星 时 钟 钟 跳 等 引起 的 数据 跳 变 异常 ， 其 中 数据 缺失 点 为 30 个 ， 如 图 
5。 最 后 以 392 个 数据 为 间隔 将 6 284 个 数据 点 分 为 16 段 ， 计 算 在 不 同 噪声 条 件 下 每 段 数据 的 多 项 式 


拟 合 残 差 标 准 差 ， 限 于 篇 幅 ， 


图 6 仅 绘 出 了 高 斯 白 噪声 标准 差 为 0. 050 时 的 拟 合 残 差 标准 差 ， 表 2 统 


计 了 两 种 观测 异常 情况 下 的 拟 合 残 差 标 准 差 与 其 均值 之 差 。 


区 段 曲线 拟 合 


— 原始 曲线 
— 多 项 式 拟 合 


弧度 
(a) 


苑 段 曲 线 拟 合 


一 一 原始 曲线 
多 项 式 拟 合 


弧度 
(b) 


图 4 多 项 式 拟 合 曲 线 与 原始 曲线 对 比 


Fig.4 Comparison between the polynomial fitting and original time-series 


4 上 数据 跳 变 
"n 


“| BKD 
1 emi, COURIR 


图 5 AMAHA RA (Er RP AEA 0. 050) 
Fig. 5 Simulated observational anomaly with White 
Gaussian Noise ( WGN) (0-0. 050) 


从 图 6 和 表 2 可 以 发 现 ， 观 测 数据 产生 缺失 
和 跳 变异 常 现 象 时 ， 多 项 式 拟 合 残 差 标准 差 较 
大 ， 这 说 明 当 载波 相位 或 伪 距 出 现 观测 异常 时 ， 
基于 多 项 式 拟 合 的 信号 评估 方法 失效 ， 因 此 必须 
截取 正常 的 观测 数据 进行 信和 号 评估 。 


4 工程 实例 


根据 信号 稳定 性 评估 方法 ， 对 北斗 系统 G02 


-e- 拟 合 残 差 标准 差 
-e- 拟 合 残 差 标准 差 均 值 


EA 
E 
> 
N 
=] 
T 


SE 


TREE 


E 6 高 斯 白 骂 声 标 准 差 为 0.050 Rr 425 FE EX 
Fig.6 Standard deviation of fitting residuals with White 
Gaussian Noise (WGN) (o=0.050) 


表 2 观测 异常 情况 下 多 项 式 拟 合 结果 
Table 2 Statistics of polynomial fitting in presence 


of observational anomaly 


高 斯 白 噪 声 。 拟 合 残 差 差 值 差 值 
标准 差 ”标准 差 均 值 (数据 缺失 ) (数据 跳 变 ) 
0. 020 0. 048 0. 176 0. 205 
0. 040 0. 065 0. 161 0. 188 
0. 050 0. 074 0. 154 0. 184 
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地 球 同 步 轨道 卫星 B2 频 点 1 支 路 2015 年 3 月 6 日 1 天 的 伪 距 测量 数据 进行 稳定 性 分 析 评 估 ， 监 测 接 
收 机 输出 数据 的 采样 间隔 为 30s， 连 续 记 录 2 880 个 观测 值 。 

图 7(a) 和 图 7(b) 分 别 绘 出 了 00: 00~24: 00 24 h 原始 伪 距 观测 值 曲线 及 其 多 项 式 拟 合 残 差 曲线 。 
由 图 7(a) 可 以 看 出 ， 该 地 球 同步 轨道 卫星 的 伪 距 观测 曲线 平滑 ， 没 有 发 生 数据 缺失 和 数据 跳 变 等 异 
常 。 从 图 7(b) 可 以 发 现 ， 伪 距 多 项 式 拟 合 残 差 曲线 在 0 附近 波动 ， 并 且 在 一 定 范围 内 变化 。 由 图 7 
(b) 还 可 以 发 现 ， 相 对 而 言 ，06: 00~09: 00 时 段 的 多 项 式 拟 合 残 差 曲 线 不 如 其 它 时 段 的 拟 合 残 差 曲 
线 变 化 平稳 。 表 3 列 出 了 以 3 h 为 一 个 时 段 划分 的 1 天 的 伪 距 拟 合 残 差 标准 差 及 其 与 标准 差 均 值 之 
差 。 从 表 3 可 以 看 出 ， 除 06: 00~09: 00 时 段 外 ， 其 它 时 段 的 伪 距 拟 合 残 差 标准 差 相 对 于 标准 差 均 值 之 
差 的 绝对 值 均 小 于 0. 100 m, 06; 00-09; 00 时 上段 的 标准 差 最 大 ， 为 0.434 m， 与 标准 差 均 值 的 偏差 为 
0. 183 m， 与 标准 差 均 值 差 值 最 小 的 时 段 为 12: 00~15: 00， 其 标准 差 为 0.218 m, 偏差 值 为 -0. 033 m, 
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È 37845 000 RE 
-0.5 
37 840 000 -10 
37 835 000 -1.5 
0 5 10 15 20 25 0 5 10 15 20 25 
时 间 / 时 间作 h 
(a) (b) 
图 7 伪 距 原始 观测 值 及 其 多 项 式 拟 合 残 差 
Fig.7 Original pseudo-range measurements and its polynomial fitting residuals 
06: 00~09: 00 时 有 段 的 伪 距 拟 合 残 差 相对 于 R3 伪 距 多 项 式 拟 合 残 差 标准 差 统 计 结 果 
其 它 时 有 段 残 差 较 大 的 原因 在 于 , 06. 00~ 09. 00 Table 3 Statistics of polynomial fitting residuals 
时 段 的 伪 距 观测 值 正好 位 于 24 h 伪 距 观测 曲线 时 要 拟 合 残 差 标 ” ”与 标准 差 均 
的 拐点 处 。 曲 线 拐点 位 置 的 拟 合 效果 一 般 没有 其 D 准 差 /m 值 之 差 /m 
它 时 有 段 的 拟 合 效 果 好 ,但 是 曲线 拐点 处 的 拟 合 不 00; 00~03: 00 0. 193 -0. 058 
会 给 整体 数据 的 拟 合 结果 带 来 影响 ， 故 多 项 式 拟 03; 00~06: 00 0. 299 0. 048 
合 结果 是 可 信和 可 用 的 ， 这 从 表 3 最 后 一 行 给 出 06: 00~09: 00 0. 434 0. 183 
的 标准 差 均 值 也 可 以 看 出 ，06: 00-09. 00 时 段 09. 00~12. 00 0. 298 0. 047 
的 拟 合 没有 影响 整体 数据 的 拟 合 结果 ， 并 且 与 12: 00~15: 00 0.218 -0. 033 
3.2 节 观 测 异 常 的 模拟 结果 相 比 ， 该 时 有 段 的 伪 距 15. 00-18. 00 0. 203 —0. 048 
拟 合 残 差 标准 差 要 比 伪 距 发 生 观 测 异 种 时 的 拟 18: 00-21. 00 0. 154 —0. 097 
合 残 差 标准 差 小 很 多 。 因而 可 以 认为 ， 06: 00~ 21: 00-24. 00 0. 212 —0. 039 
09: 00 时 段 的 拟 合 残 差 表现 较 大 的 现象 与 拟 合 时 均值 0. 251 0. 000 


段 划 分 有 关 ， 而 不 是 由 信号 稳定 性 起 伏 引 起 的 。 

为 了 进一步 分 析 06: 00~ 09 : 00 时 段 的 伪 距 拟 合 残 差 结果 ， 对 残 差 统计 特性 进行 讨论 。 多 项 式 拟 
合 残 差 主要 包含 混杂 在 原始 数据 中 的 随机 噪声 成 分 ， 换 言 之 ， 伪 距 拟 合 残 差 通常 是 由 观测 噪声 产生 的 
随机 波动 ， 理 论 上 应 该 服从 零 均 值 的 高 斯 分 布 ， 然 而 由 (1) 式 可 知 ， 卫 星 信号 在 空间 传播 过 程 中 受到 
对 流 层 和 电离 层 折射 误差 等 多 种 因素 的 影响 ， 从 而 导致 总 体 方差 和 总 体 均 值 产 生变 化 ， 若 想 知 道 总 体 
方差 和 总 体 均值 的 变化 情况 ， 需 要 对 两 个 高 斯 分 布 总 体 方差 比 或 总 体 均 值 差 进行 估计 。 这 里 仅 估计 两 
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个 高 斯 分 布 总 体 均 值 差 。 
选取 06: 00~09: 00 和 15: 00~18: 00 两 个 时 段 的 多 项 式 拟 合 残 差 值 (概率 分 布 如 图 8) 作为 两 
个 总 体 ( 方 差 和 均值 等 参数 见 表 4)， 估计 两 个 总 体 均值 差 的 一 个 置信 和 度 为 1-a 的 置信 和 区间， 


[on [on We x . 
0,7 -HtZ /一 + 一 ， 其 中 ，Z 表示 标准 高 斯 分 布 的 上 o/2 分 位 点 ， 本 文 取 w=0.05。 经 计算 


1 ni 
获得 两 个 总 体 均 值 差 的 95% 置 信 区 间 为 [ -0. 358，0. 237] 。 因 为 此 置信 区 间 包 括 0， 因 而 可 以 判断 这 
两 个 时 段 伪 距 拟 合 残 差 的 均值 没有 显著 差别 。 


à 


0.9 1.0 
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0.6 0.7 
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E 04 *X 
Sue E 04 
im 0.3 
0.2 02 
0.1 0.1 
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-2 -1 0 1 2 -1 -0.5 0 0.5 1 
拟 合 残 差 幅 度 /m 拟 合 残 差 幅度 /m 
(a) 06:00~09:00 (b) 15:00-18:00 
图 8 伪 距 多 项 式 拟 合 残 差 概率 分 布 
Fig. 8 Probability distribution of polynomial fitting residuals 
5 结束 语 RA ， 伪 距 多 项 式 拟 合 残 差 统计 参数 
Table 4 Statistic parameters of polynomial fitting 
为 了 分 析 卫星 导航 信号 的 稳定 性 ， 以 最 小 二 pep 


乘 拟 合 原理 为 基础 ， 通 过 对 监测 接收 机 的 载波 相 高 斯 总 体 样本 数 总 体 均值 。 总 体 方差 
位 和 伪 距 测量 数据 进行 评估 ， 实 现 了 对 全 球 卫 星 06:00-09:00  n,210 j=0.040  g-0.189 
导航 系统 空间 信号 稳定 性 的 分 析 。 数 字 仿 真 结 15:00~18:00 m,=10 p=0.021  g;-0.041 
表明 ， 多 项 式 拟 合 方法 在 曲线 拐点 处 的 拟 合 效 果 ” 
不 如 曲线 其 它 处 ,但 曲线 拐点 处 的 拟 合 不 会 影响 整体 数据 的 拟 合 结果 。 当 载波 相位 或 伪 距 观测 量 因 接 
收 机 失 锁 等 因素 出 现 数据 缺失 或 因 卫 星 时 钟 等 因素 引起 数据 跳 变 观 测 异 常 时 ， 多 项 式 拟 合法 不 再 适 
用 ， 因 而 必须 选取 正常 观测 数据 进行 分 析 。 通 过 对 北斗 系统 地 球 同步 轨道 卫星 伪 距 实测 数据 的 评估 可 
以 发 现 ， 该 卫星 发 射 的 I 支 路 信号 在 1 天 内 是 稳定 的 。 由 于 多 项 式 拟 合 方法 具有 一 定 的 局 限 性 ， 给 信 
号 评估 结果 带 来 一 定 的 影响 ， 因 此 如 何 对 拟 合 方法 进行 改进 是 下 一 步 的 研究 重点 。 此 外 ， 如 何 根 据 卫 
星 信 号 的 变化 趋势 进行 信号 的 长 期 稳定 性 评估 也 是 未 来 研究 的 关注 点 。 
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Study on Methods for Analyzing the Stability of Satellite Navigation Signals 
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Abstract: The quality assessment of GNSS ( Global Navigation Satellite System) space signals plays an 


important role in the design and operation of the system. It is also an essential task in the integrity monitoring 


of space signals. In this paper, the methodology for evaluating GNSS space signals is first described in detail, 


which is based on pseudo-range or carrier-phase observations from a GNSS monitoring receiver. Next, the 


feasibility of the presented algorithm is validated by using a simulation experiment. Furthermore, an 


experiment on the real BDS pseudo-range measurements is carried out in order to validate the effectiveness of 


the developed method. The results show that the proposed scheme is highly potential to evaluate the quality of 


GNSS space signals according to the three metrics, namely original measurements, fitting residuals and 


standard deviation of fitting residuals. Therefore, the proposed evaluation technique can provide a guarantee 


and help for the continuity and reliability of GNSS space signals. 


Key words: Global Navigation Satellite System (GNSS) ; Analysis of the signal quality; Stability; Pseudo-range 


